Rapid genome wide mapping of phosphine resistance loci by a simple regional averaging analysis in the red flour beetle, Tribolium castaneum

Background Next-generation sequencing technology is an important tool for the rapid, genome-wide identification of genetic variations. However, it is difficult to resolve the ‘signal’ of variations of interest and the ‘noise’ of stochastic sequencing and bioinformatic errors in the large datasets that are generated. We report a simple approach to identify regional linkage to a trait that requires only two pools of DNA to be sequenced from progeny of a defined genetic cross (i.e. bulk segregant analysis) at low coverage (<10×) and without parentage assignment of individual SNPs. The analysis relies on regional averaging of pooled SNP frequencies to rapidly scan polymorphisms across the genome for differential regional homozygosity, which is then displayed graphically. Results Progeny from defined genetic crosses of Tribolium castaneum (F4 and F19) segregating for the phosphine resistance trait were exposed to phosphine to select for the resistance trait while the remainders were left unexposed. Next generation sequencing was then carried out on the genomic DNA from each pool of selected and unselected insects from each generation. The reads were mapped against the annotated T. castaneum genome from NCBI (v3.0) and analysed for SNP variations. Since it is difficult to accurately call individual SNP frequencies when the depth of sequence coverage is low, variant frequencies were averaged across larger regions. Results from regional SNP frequency averaging identified two loci, tc_rph1 on chromosome 8 and tc_rph2 on chromosome 9, which together are responsible for high level resistance. Identification of the two loci was possible with only 5-7× average coverage of the genome per dataset. These loci were subsequently confirmed by direct SNP marker analysis and fine-scale mapping. Individually, homozygosity of tc_rph1 or tc_rph2 results in only weak resistance to phosphine (estimated at up to 1.5-2.5× and 3-5× respectively), whereas in combination they interact synergistically to provide a high-level resistance >200×. The tc_rph2 resistance allele resulted in a significant fitness cost relative to the wild type allele in unselected beetles over eighteen generations. Conclusion We have validated the technique of linkage mapping by low-coverage sequencing of progeny from a simple genetic cross. The approach relied on regional averaging of SNP frequencies and was used to successfully identify candidate gene loci for phosphine resistance in T. castaneum. This is a relatively simple and rapid approach to identifying genomic regions associated with traits in defined genetic crosses that does not require any specialised statistical analysis.

Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia stem rot in Brassica napus

Brassica napus is one of the most important oil crops in the world, and stem rot caused by the fungus Sclerotinia sclerotiorum results in major losses in yield and quality. To elucidate resistance genes and pathogenesis-related genes, genome-wide association analysis of 347 accessions was performed using the Illumina 60K Brassica SNP (single nucleotide polymorphism) array. In addition, the detached stem inoculation assay was used to select five highly resistant (R) and susceptible (S) B. napus lines, 48 h postinoculation with S. sclerotiorum for transcriptome sequencing. We identified 17 significant associations for stem resistance on chromosomes A8 and C6, five of which were on A8 and 12 on C6. The SNPs identified on A8 were located in a 409-kb haplotype block, and those on C6 were consistent with previous QTL mapping efforts. Transcriptome analysis suggested that S. sclerotiorum infection activates the immune system, sulphur metabolism, especially glutathione (GSH) and glucosinolates in both R and S genotypes. Genes found to be specific to the R genotype related to the jasmonic acid pathway, lignin biosynthesis, defence response, signal transduction and encoding transcription factors. Twenty-four genes were identified in both the SNP-trait association and transcriptome sequencing analyses, including a tau class glutathione S-transferase (GSTU) gene cluster. This study provides useful insight into the molecular mechanisms underlying the plant's response to S. sclerotiorum.

Stem cankers on sunflower (Helianthus annuus) in Australia reveal a complex of pathogenic Diaporthe (Phomopsis) species.

The identification of Diaporthe (anamorph Phomopsis) species associated with stem canker of sunflower (Helianthus annuus) in Australia was studied using morphology, DNA sequence analysis and pathology. Phylogenetic analysis revealed three clades that did not correspond with known taxa, and these are believed to represent novel species. Diaporthe gulyae sp. nov. is described for isolates that caused a severe stem canker, specifically pale brown to dark brown, irregularly shaped lesions centred at the stem nodes with pith deterioration and mid-stem lodging. This pathogenicity of D. gulyae was confirmed by satisfying Koch's Postulates. These symptoms are almost identical to those of sunflower stem canker caused by D. helianthi that can cause yield reductions of up to 40% in Europe and the USA, although it has not been found in Australia. We show that there has been broad misapplication of the name D. helianthi to many isolates of Diaporthe ( Phomopsis) found causing, or associated with, stem cankers on sunflower. In GenBank, a number of isolates had been identified as D. helianthi, which were accommodated in several clades by molecular phylogenetic analysis. Two less damaging species, D. kochmanii sp. nov. and D. kongii sp. nov., are also described from cankers on sunflower in Australia.

Spot form of net blotch resistance in barley is under complex genetic control

Key message: Evaluation of resistance toPyrenophora teresf.maculatain barley breeding populations via association mapping revealed a complex genetic architecture comprising a mixture of major and minor effect genes. Abstract: In the search for stable resistance to spot form of net blotch (Pyrenophora teres f. maculata, SFNB), association mapping was conducted on four independent barley (Hordeum vulgare L.) breeding populations comprising a total of 898 unique elite breeding lines from the Northern Region Barley Breeding Program in Australia for discovery of quantitative trait loci (QTL) influencing resistance at seedling and adult plant growth stages. A total of 29 significant QTL were validated across multiple breeding populations, with 22 conferring resistance at both seedling and adult plant growth stages. The remaining 7 QTL conferred resistance at either seedling (2 QTL) or adult plant (5 QTL) growth stages only. These 29 QTL represented 24 unique genomic regions, of which five were found to co-locate with previously identified QTL for SFNB. The results indicated that SFNB resistance is controlled by a large number of QTL varying in effect size with large effects QTL on chromosome 7H. A large proportion of the QTL acted in the same direction for both seedling and adult responses, suggesting that phenotypic selection for SFNB resistance performed at either growth stage could achieve adequate levels of resistance. However, the accumulation of specific resistance alleles on several chromosomes must be considered in molecular breeding selection strategies.

Molecular Breeding for Complex Adaptive Traits: How Integrating Crop Ecophysiology and Modelling Can Enhance Efficiency

Progress in crop improvement is limited by the ability to identify favourable combinations of genotypes (G) and management practices (M) in relevant target environments (E) given the resources available to search among the myriad of possible combinations. To underpin yield advance we require prediction of phenotype based on genotype. In plant breeding, traditional phenotypic selection methods have involved measuring phenotypic performance of large segregating populations in multi-environment trials and applying rigorous statistical procedures based on quantitative genetic theory to identify superior individuals. Recent developments in the ability to inexpensively and densely map/sequence genomes have facilitated a shift from the level of the individual (genotype) to the level of the genomic region. Molecular breeding strategies using genome wide prediction and genomic selection approaches have developed rapidly. However, their applicability to complex traits remains constrained by gene-gene and gene-environment interactions, which restrict the predictive power of associations of genomic regions with phenotypic responses. Here it is argued that crop ecophysiology and functional whole plant modelling can provide an effective link between molecular and organism scales and enhance molecular breeding by adding value to genetic prediction approaches. A physiological framework that facilitates dissection and modelling of complex traits can inform phenotyping methods for marker/gene detection and underpin prediction of likely phenotypic consequences of trait and genetic variation in target environments. This approach holds considerable promise for more effectively linking genotype to phenotype for complex adaptive traits. Specific examples focused on drought adaptation are presented to highlight the concepts.

Refined Global Analysis of Bemisia tabaci (Hemiptera: Sternorrhyncha: Aleyrodoidea: Aleyrodidae) Mitochondrial Cytochrome Oxidase 1 to Identify Species Level Genetic Boundaries.

Identifying species boundaries within morphologically indistinguishable cryptic species complexes is often contentious. For the whitefly Bemisia tabaci (Gennadius) (Hemiptera: Sternorrhyncha: Aleyrodoidea: Aleyrodidae), the lack of a clear understanding about the genetic limits of the numerous genetic groups and biotypes so far identified has resulted in a lack of consistency in the application of the terms, the approaches use to apply them and in our understanding of what genetic structure within B. tabaci means. Our response has been to use mitochondrial gene cytochrome oxidase one to consider how to clearly and consistently define genetic separation. Using Bayesian phylogenetic analysis and analysis of sequence pairwise divergence we found a considerably higher to number of genetic groups than had been previously determined with two breaks in the distribution, one at 11% and another at 3.5%. At >11% divergence, 11 distinct groups were resolved, whereas at >3.5% divergence 24 groups were identified. Consensus sequences for each of these groups were determined and were shown to be useful in the correct assignment of sequences of unknown origin. The 3.5% divergence bound is consistent with species level separations in other insect taxa and Suggests that B. tabaci is it cryptic species composed of at least 24 distinct species. We further show that the placement of Bemesia atriplex (Froggatt) within the B. tabaci in, group adds further weight to the argument for species level separation within B. tabaci. This new analysis, which constructs consensus sequences and uses these its a standard against which unknown sequences call be compared, provides for the first time it consistent means of identifying the genetic hounds of each species with it high degree of certainty.

Gene to Phenotype Integrative Science and Modelling Applied to Pioneer elite

This project aims to use simulatiion modelling to improve our understanding of the genetics and physiology of complex traits with a view to increasing the rate of genetic gain in plant breeding programs.

Evidence for different QTL underlying the immune and hypersensitive responses of Eucalyptus globulus to the rust pathogen Puccinia psidii

The rust Puccinia psidii infects many species in the family Myrtaceae. Native to South America, the pathogen has recently entered Australia which has a rich Myrtaceous flora, including trees of the ecologically and economically important genus Eucalyptus. We studied the genetic basis of variation in rust resistance in Eucalyptus globulus, the main plantation eucalypt in Australia. Quantitative trait loci (QTL) analysis was undertaken using 218 genotypes of an outcross F2 mapping family, phenotyped by controlled inoculation of their open pollinated progeny with the strain of P. psidii found in Australia. QTL analyses were conducted using a binary classification of individuals with no symptoms (immune) versus those with disease symptoms, and in a separate analysis dividing plants with disease symptoms into those exhibiting the hypersensitive response versus those with more severe symptoms. Four QTL were identified, two influencing whether a plant exhibited symptoms (Ppr2 and Ppr3), and two influencing the presence or absence of a hypersensitive reaction (Ppr4 and Ppr5). These QTL mapped to four different linkage groups, none of which overlap with Ppr1, the major QTL previously identified for rust resistance in Eucalyptus grandis. Candidate genes within the QTL regions are presented and possible mechanisms discussed. Together with past findings, our results suggest that P. psidii resistance in eucalypts is quantitative in nature and influenced by the complex interaction of multiple loci of variable effect.

Domestication to crop improvement: Genetic resources for Sorghum and Saccharum (Andropogoneae).

Background: Both sorghum (Sorghum bicolor) and sugarcane (Saccharum officinarum) are members of the Andropogoneae tribe in the Poaceae and are each other's closest relatives amongst cultivated plants. Both are relatively recent domesticates and comparatively little of the genetic potential of these taxa and their wild relatives has been captured by breeding programmes to date. This review assesses the genetic gains made by plant breeders since domestication and the progress in the characterization of genetic resources and their utilization in crop improvement for these two related species. Genetic Resources: The genome of sorghum has recently been sequenced providing a great boost to our knowledge of the evolution of grass genomes and the wealth of diversity within S. bicolor taxa. Molecular analysis of the Sorghum genus has identified close relatives of S. bicolor with novel traits, endosperm structure and composition that may be used to expand the cultivated gene pool. Mutant populations (including TILLING populations) provide a useful addition to genetic resources for this species. Sugarcane is a complex polyploid with a large and variable number of copies of each gene. The wild relatives of sugarcane represent a reservoir of genetic diversity for use in sugarcane improvement. Techniques for quantitative molecular analysis of gene or allele copy number in this genetically complex crop have been developed. SNP discovery and mapping in sugarcane has been advanced by the development of high-throughput techniques for ecoTILLING in sugarcane. Genetic linkage maps of the sugarcane genome are being improved for use in breeding selection. The improvement of both sorghum and sugarcane will be accelerated by the incorporation of more diverse germplasm into the domesticated gene pools using molecular tools and the improved knowledge of these genomes.

SSR analysis of cultivated groundnut (Arachis hypogaea L.) germplasm resistant to rust and late leaf spot diseases

Cultivated groundnut (Arachis hypogaea L.) is an agronomically and economically important oilseed crop grown extensively throughout the semi-arid tropics of Asia, Africa and Latin America. Rust (Puccinia arachidis) and late leaf spot (LLS, Phaseoisariopsis personata) are among the major diseases causing significant yield loss in groundnut. The development of varieties with high levels of resistance has been constrained by adaptation of disease isolates to resistance sources and incomplete resistance in resistant sources. Despite the wide range of morphological diversity observed in the cultivated groundnut gene pool, molecular marker analyses have thus far been unable to detect a parallel level of genetic diversity. However, the recent development of simple sequence repeat (SSR) markers presents new opportunities for molecular diversity analysis of cultivate groundnut. The current study was conducted to identify diverse disease resistant germplasm for the development of mapping populations and for their introduction into breeding programs. Twenty-three SSRs were screened across 22 groundnut genotypes with differing levels of resistance to rust and LLS. Overall, 135 alleles across 23 loci were observed in the 22 genotypes screened. Twelve of the 23 SSRs (52%) showed a high level of polymorphism, with PIC values ≥0.5. This is the first report detecting such high levels of genetic polymorphism in cultivated groundnut. Multi-dimensional scaling and cluster analyses revealed three well-separated groups of genotypes. Locus by locus AMOVA and Kruskal-Wallis one-way ANOVA identified candidate SSR loci that may be valuable for mapping rust and LLS resistance. The molecular diversity analysis presented here provides valuable information for groundnut breeders designing strategies for incorporating and pyramiding rust and late leaf spot resistances and for molecular biologists wishing to create recombinant inbred line populations to map these traits.

Identification of diverse full-length endogenous betaretroviruses in megabats and microbats

Background: Betaretroviruses infect a wide range of species including primates, rodents, ruminants, and marsupials. They exist in both endogenous and exogenous forms and are implicated in animal diseases such as lung cancer in sheep, and in human disease, with members of the human endogenous retrovirus-K (HERV-K) group of endogenous betaretroviruses (βERVs) associated with human cancers and autoimmune diseases. To improve our understanding of betaretroviruses in an evolutionarily distinct host species, we characterized βERVs present in the genomes and transcriptomes of mega- and microbats, which are an important reservoir of emerging viruses.Results: A diverse range of full-length βERVs were discovered in mega- and microbat genomes and transcriptomes including the first identified intact endogenous retrovirus in a bat. Our analysis revealed that the genus Betaretrovirus can be divided into eight distinct sub-groups with evidence of cross-species transmission. Betaretroviruses are revealed to be a complex retrovirus group, within which one sub-group has evolved from complex to simple genomic organization through the acquisition of an env gene from the genus Gammaretrovirus. Molecular dating suggests that bats have contended with betaretroviral infections for over 30 million years.Conclusions: Our study reveals that a diverse range of betaretroviruses have circulated in bats for most of their evolutionary history, and cluster with extant betaretroviruses of divergent mammalian lineages suggesting that their distribution may be largely unrestricted by host species barriers. The presence of βERVs with the ability to transcribe active viral elements in a major animal reservoir for viral pathogens has potential implications for public health. © 2013 Hayward et al.; licensee BioMed Central Ltd.

An inordinate fondness for Fusarium: Phylogenetic diversity of fusaria cultivated by ambrosia beetles in the genus Euwallacea on avocado and other plant hosts

Ambrosia beetle fungiculture represents one of the most ecologically and evolutionarily successful symbioses, as evidenced by the 11 independent origins and 3500 species of ambrosia beetles. Here we document the evolution of a clade within Fusarium associated with ambrosia beetles in the genus Euwallacea (Coleoptera: Scolytinae). Ambrosia Fusarium Clade (AFC) symbionts are unusual in that some are plant pathogens that cause significant damage in naive natural and cultivated ecosystems, and currently threaten avocado production in the United States, Israel and Australia. Most AFC fusaria produce unusual clavate macroconidia that serve as a putative food source for their insect mutualists. AFC symbionts were abundant in the heads of four Euwallacea spp., which suggests that they are transported within and from the natal gallery in mandibular mycangia. In a four-locus phylogenetic analysis, the AFC was resolved in a strongly supported monophyletic group within the previously described Cade 3 of the Fusarium solani species complex (FSSC). Divergence-time estimates place the origin of the AFC in the early Miocene similar to 21.2 Mya, which coincides with the hypothesized adaptive radiation of the Xyleborini. Two strongly supported clades within the AFC (Clades A and B) were identified that include nine species lineages associated with ambrosia beetles, eight with Euwallacea spp. and one reportedly with Xyleborus ferrugineus, and two lineages with no known beetle association. More derived lineages within the AFC showed fixation of the clavate (club-shaped) macroconidial trait, while basal lineages showed a mix of clavate and more typical fusiform macroconidia. AFC lineages consisted mostly of genetically identical individuals associated with specific insect hosts in defined geographic locations, with at least three interspecific hybridization events inferred based on discordant placement in individual gene genealogies and detection of recombinant loci. Overall, these data are consistent with a strong evolutionary trend toward obligate symbiosis coupled with secondary contact and interspecific hybridization. (C) 2013 Elsevier Inc. All rights reserved.

Rapid assessment of genotype-by-environment interactions and heritability for growth rate in aquaculture species using in vitro fertilisation and DNA tagging

Commercial environments may receive only a fraction of expected genetic gains for growth rate as predicted from the selection environment. This fraction is result of undesirable genotype-by-environment interactions (GxE) and measured by the genetic correlation (rg) of growth between environments. Rapid estimates of genetic correlation achieved in one generation are notoriously difficult to estimate with precision. A new design is proposed where genetic correlations can be estimated by utilising artificial mating from cryopreserved semen and unfertilised eggs stripped from a single female. We compare a traditional phenotype analysis of growth to a threshold model where only the largest fish are genotyped for sire identification. The threshold model was robust to differences in family mortality differing up to 30%. The design is unique as it negates potential re-ranking of families caused by an interaction between common maternal environmental effects and growing environment. The design is suitable for rapid assessment of GxE over one generation with a true 0.70 genetic correlation yielding standard errors as low as 0.07. Different design scenarios were tested for bias and accuracy with a range of heritability values, number of half-sib families created, number of progeny within each full-sib family, number of fish genotyped, number of fish stocked, differing family survival rates and at various simulated genetic correlation levels.

Phylodynamic evidence of the migration of turnip mosaic potyvirus from Europe to Australia and New Zealand

Turnip mosaic virus (TuMV) is a potyvirus that is transmitted by aphids and infects a wide range of plant species. We investigated the evolution of this pathogen by collecting 32 isolates of TuMV, mostly from Brassicaceae plants, in Australia and New Zealand. We performed a variety of sequence-based phylogenetic and population genetic analyses of the complete genomic sequences and of three non-recombinogenic regions of those sequences. The substitution rates, divergence times and phylogeographical patterns of the virus populations were estimated. Six inter- and seven intralineage recombination-type patterns were found in the genomes of the Australian and New Zealand isolates, and all were novel. Only one recombination-type pattern has been found in both countries. The Australian and New Zealand populations were genetically different, and were different from the European and Asian populations. Our Bayesian coalescent analyses, based on a combination of novel and published sequence data from three nonrecombinogenic protein-encoding regions, showed that TuMV probably started to migrate from Europe to Australia and New Zealand more than 80 years ago, and that distinct populations arose as a result of evolutionary drivers such as recombination. The basal-B2 subpopulation in Australia and New Zealand seems to be older than those of the world-B2 and -B3 populations. To our knowledge, our study presents the first population genetic analysis of TuMV in Australia and New Zealand. We have shown that the time of migration of TuMV correlates well with the establishment of agriculture and migration of Europeans to these countries.